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We propose a new mechanism of leptogenesis in which the asymmetries in lepton numbers are 
produced through the CP-violating oscillations of "sterile" (electroweak singlet) neutrinos. The 
asymmetry is communicated from singlet neutrinos to ordinary leptons through their Yukawa cou- 
plings. The lepton asymmetry is then reprocessed into baryon asymmetry by electroweak sphalerons. 
We show that the observed value of baryon asymmetry can be generated in this way, and the masses 
of ordinary neutrinos induced by the seesaw mechanism are in the astrophysically and cosmologically 
interesting range. Except for singlet neutrinos, no physics beyond the Standard Model is required. 
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1. The origin of the excess of baryons over anti-baryons 
in the Universe remains one of the fascinating problems of 
particle physics and cosmology. A number of mechanisms 
have been proposed to date to explain this asymmetry 
(for recent reviews see, e.g., 0). One of the simplest pos- 
sibilities, suggested by Fukugita and Yanagida is that 
the baryon asymmetry has originated from physics in the 
leptonic sector. Namely, it was assumed that at temper- 
atures well above the electroweak scale, lepton asymme- 
try was produced, which was then reprocessed into the 
baryon asymmetry by non-perturbative electroweak ef- 
fects [|| - sphalerons Q . According to ref . || the lepton 
asymmetry is generated in out-of-equilibrium, CP- and 
lepton number non-conserving decays of heavy Majorana 
neutrinos (for recent discussions see, e.g., ref. Q and ref- 
erences therein). 

In this Letter we propose a new realization of baryo- 
genesis through leptogenesis which also makes use of 
the electroweak reprocessing of the lepton number into 
the baryon number. Like the Fukugita- Yanagida mech- 
anism, our proposal requires only mild extension of 
the Standard Model by introducing "sterile" (i.e., elec- 
troweak singlet) heavy neutrinos. However, our mecha- 
nism of leptogenesis is entirely different from that of ref. 
Q : we suggest that asymmetries in lepton numbers were 
generated due to oscillations of these singlet neutrinos 
and their interactions with ordinary matter in the early 
Universe. Moreover, the novel feature of our scenario 
is that the total lepton number is not violated in these 
oscillations and/or interactions; an important ingredient 
is separation (rather than generation) of lepton number, 
i.e., its redistribution between different species of singlet 
neutrinos. 

For this reason we do not necessarily require that sin- 
glet neutrinos be Majorana particles; Dirac "sterile" neu- 
trinos are equally suitable (and even better in some re- 
spect) for our mechanism to work. Furthermore, in our 
case the values of the masses and couplings of singlet 



neutrinos are very different from those of ref. [|2j. 

2. Let us consider the Standard Model extended 
by adding three types of Majorana neutrinos N a , a — 
A, B, C which interact with other particles only through 
their Yukawa couplings The corresponding La- 

grangian can be written in the "Yukawa basis" (where 
the matrix of Yukawa coupling constants has been diag- 
onalized) as follows, 



C = N Ra i@N Ra + h a l a N Ra 3> 



Mr, 



N Ra CN Rb 



Here N Ra are right-handed components of N a , h a are 
(real) Yukawa coupling constants, l a are three indepen- 
dent linear combinations of usual leptonic doublets, $ is 
the Higgs doublet, and M a b is the mass matrix. 

We suggest that the baryogenesis proceeds in the fol- 
lowing way: 

(i) In the course of the evolution of the Universe, sin- 
glet neutrinos are produced through their Yukawa cou- 
plings. The production mechanism of singlet neutrinos 
conserves CP, i.e. for each type equal numbers of par- 
ticles and antiparticles (particles of opposite helicities in 
the Majorana case) are produced. 

(ii) Once created, singlet neutrinos oscillate, and also 
interact with ordinary matter. None of these processes 
violates the total lepton number L tot = L+La+Lb+Lq, 
where L is the usual lepton number (we assume that Ma- 
jorana masses are small enough, see below). However, CP 
is not conserved due to mixing in the singlet neutrino sec- 
tor. Therefore the initially created state with individual 
lepton numbers La = L R = Lq = evolves through the 
oscillations into a state in which La ^ 0, L R ^ 0, Lq ^ 
but still L tot = 0. That is, the total lepton number gets 
unevenly distributed between different species. 

(iii) Singlet neutrinos communicate their lepton asym- 
metry to ordinary neutrinos and charged leptons through 
their Yukawa couplings. We assume that the Yukawa 
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couplings of Na, Nb and Nc have certain hierarchy, so 
that neutrinos of at least one type, Na , come into ther- 
mal equilibrium before the time tEW at which sphalerons 
become inoperative (the corresponding temperature is 
Tew ~ 100 GeV) and those of at least one other type, 
Nc, do not equilibrate by tEW ■ The neutrinos of the 
third type, Nb, may or may not come into thermal equi- 
librium by tEw- To be specific, we discuss mostly the 
case when the Yukawa coupling of Nb is relatively large 
so that Nb equilibrate at temperatures exceeding Tew 
(the opposite case is treated in a similar way). In this 
case the lepton numbers La and Lb are communicated 
to the ordinary leptons before tEw whereas Lc is not. 
Therefore (a fraction of) La + Lb is reprocessed into 
baryon asymmetry by electroweak sphalerons; the lepton 
number Lc is transfered to the active leptons only af- 
ter the sphalerons have already switched off, so it has no 
effect on the baryogenesis. 

Let us stress that if not for the sphaleron freezing ef- 
fect, no baryon asymmetry would have been obtained: all 
three singlet neutrino species eventually transform into 
ordinary leptons, and since L tot = 0, no net lepton charge 
would have been generated in the sector of ordinary lep- 
tons. Hence, the requirement that Nc do not get into 
thermal equilibrium before tEW is crucial for our mech- 
anism. 

3. The system of three types of singlet neutrino of a 
given momentum k{t) oc T(t), interacting with cosmic 
plasma, is described by 3 x 3 density matrix pit) which 
obeys the evolution equation [Q 



dp 
i — 
dt 



[h,p]-Ut, p } + -{tp,i-p} 



(i) 



Here H is Hermitean effective Hamiltonian, V and T p are 
destruction and production rates, respectively, and we 
have neglected the relaxation due to elastic scattering, 
see below. In what follows we will use the approximation 
of Boltzmann statistics for order-of-magnitude analysis; 
in this approximation the last term in eq. (jl|) becomes 
merely iT p . 

In the Yukawa basis one has at T ^> M a 



M 2 



(2) 



Here U is the mixing matrix which relates the Yukawa 
states and the mass eigenstates Ni [i = 1,2,3): 
N a = U a iNi (we assume that mixing is small, and 
choose the standard parametrization for U §); M 2 = 
diag {Ml, Ml, Ml) is the matrix of mass eigenvalues; 
k(t) is the neutrino momentum which depends on time 
due to the expansion of the Universe. The CP violation 
in the system is described by the CP-odd phase 5 in the 
mixing matrix U. 

The medium effects are described by real potential V 
and rates T and T p — exp(— fc/T)r || whose matrices 
are diagonal in the Yukawa basis, 



v = dm g (v A ,v B ,Vc), r = dia g (r A ,r B ,r c ) 

The potentials V a are due to the coherent forward scatter- 
ing processes, the main contribution coming from the 1- 
loop self energy diagrams with ordinary lepton and Higgs 
doublets in the intermediate states. For our estimates we 
use their thermal averages which at temperatures above 
T EW are p| 



V a = \h\T. 



(3) 



Consider now the rates T a . We will be interested in the 
temperatures far exceeding the masses of singlet neutri- 
nos. Therefore the rates of the 1 <-> 2 reactions which 
correspond to the absorptive parts of the self energy dia- 
grams are suppressed by the factor M a /T and 2^2 reac- 
tions are more important. The main contributions to r a 
come from the Higgs exchange reactions QlNro. tnl a , 
t£N Ra <-> Q£l a and l a c N Ra <-> t R Q£ where Q L is the 
third-generation quark doublet. The result |ll]] for the 
average destruction rates at T ^> M a (corrected to in- 
clude the color factor for quarks) is 



r„ 



9hl 
64tt 3 



KT. 



(4) 



Here h t ~ 1 is the top quark Yukawa coupling. The rates 
of the elastic 2 — > 2 scattering processes are proportional 
to h 2 hl instead of h 2 h 2 . We will need very small Yukawa 
couplings of singlet neutrinos, so the elastic processes can 
be safely neglected. 

At late times the lepton asymmetry is stored in the 
least interacting species of singlet neutrino, which at 
t ~ tEw coincides with the mass eigenstate N$ w Nc- 
The conditions that iVi(~ Na) and iV 2 (« Ng) come into 
thermal equilibrium before the time tEW, while N3 do 
not, are 

r M (7W) > h(t ew ) , r 3 (iw) < h{t ew ) , (5) 

where H(T) = T 2 /M^ is the Hubble parameter, M£, = 
M P1 /1.66vk - 10 18 GeV, and T 3 sa Tc + s 2 13 r A + 4 3 TB- 
Here S13 = sin #13 and S23 = sin #23 determine the admix- 
tures of Na and Nb in V 3 . 

The conditions (||) translate into bounds on the 
Yukawa couplings, 



L A,B 



> 2 • 10" 



and on mixing angles, 

sj 3 < 2 • 1Q- I4 h- A 2 



A ) 



h%<2- 10" 14 . 



s| 3 < 2 • 10- 14 ^ 2 



(0) 



(7) 



Eq. (||) implies a certain hierarchy between the Yukawa 
couplings, which, however, need not be very strong. 

4. As described above, to find the baryon asymmetry 
we should calculate the asymmetry AL(tEw) which was 
communicated to usual leptons by the time tsw at which 
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sphalerons switch off. This asymmetry emerges because 
singlet neutrinos N$ do not transfer their asymmetry to 
active neutrinos by the time t ew due to smallness of he, 
S13 and S23. Since the total lepton number is conserved in 
all processes of interest, we have Al^ew) = ~ Aq^ew) 
(up to a factor of order one that accounts for the distribu- 
tion of the asymmetry between B, L and La+Lb', we will 
not write this factor in formulas below), where A 3 (tEw) 
is the asymmetry stored in N3. The asymmetry A 3 (tEw) 
can be found as follows. Let S(t, to) be the evolution ma- 
trix corresponding to the operator H = H — (i/2)T (no- 
tice that S(t, to) is not unitary since this operator is non- 
Hermitean). The density matrix p(t) can be expressed 
through S(t,to) assuming p(ti) — where U is the time 
at which the production of singlet neutrino begins Jl2| . 
The ratio of the number density of N3, 113 = ^33, to the 
equilibrium density of one spin degree of freedom at time 
tEw is 



n 3 (t E w) 
n eq (tEw) 



a , b 



= J2 dt ° r a{to)\U f 3b S ba {tEW,to)\ 2 , 



(8) 



The asymmetry A 3 (tEw) is the CP-odd part of the quan- 
tity®. 

The integration over the production time to can be 
performed in a closed form. Indeed, the matrices S and 
St obey d to S = iSH(t ), d to S^ = -iH^(t )S^. From 
these equations and T — i(H — H>), one finds that 
dt (SS') = ST(t )S^. Using this relation one can read- 
ily perform the integration over the production time to 
inEq. (ph, 



n 3 (t 



EW) 



' eq 



(tEw) 



— 1 — [S (tEW, U)S M (tEW, U) 
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Here S = U'SU is the evolution matrix in the mass 
eigenstate basis. The CP-odd part of this expression de- 
termines the asymmetry transfered to usual leptons by 
tEW (and hence the generated baryon asymmetry) 



A L (t E w) = (n L - n L )/n 7 = -A 3 (t E w) 



9 ^2 \ S M( t EW,U)\ 



CP-odd 1 



(9) 



where the factor 1/2 accounts for two helicity states of 
photon. The production of singlet neutrinos starts at 
very early times, so we set U = 0. Since Tew is much 
smaller than all relevant energy parameters in the prob- 
lem, one can formally let tsw — > 00 in actual calcula- 
tions. 

5. The lepton asymmetry is produced mainly at the 
epoch tL when the differences of the eigenvalues of the 
Hamiltonian fiy = Qj— Qj become of order of the Hubble 
parameter H: Oy(t/,) ~ 1/tz,- Indeed, at t <C iz, the 
elements of p(t) essentially stay constant, whereas at t ^> 



tL they undergo fast oscillations and due to averaging 
effects the asymmetry is strongly suppressed. 

In what follows we will present the results for the most 
interesting range of the parameter space where the mass 
differences are relatively large, AM 2 > (h 2 lb /8) 3 M* P ? 
[h\ h = h\- hf), AM 2 being the typical value of AM?-. 
In this case the mass terms dominate over potentials in 
H at the epoch t^, so that ~ AM?-/2T and the lep- 

togenesis temperature is T L = T(t L ) ~ (M^AM?) 1 / 3 . 
At t = ti, we have 



H 8 \ AM 2 / 



1/3 



A< 1. 



(10) 



This means that the potentials V a and rates T a can be 
treated perturbatively with A being the expansion pa- 
rameter. The lepton asymmetry ([}]) appears in the third 
order of perturbation theory, and therefore is suppressed 
by the cube of A. This can be seen in the mass eigenstate 
basis, where the potential has the form WVU. Indeed, 
Ai is a CP - violating observable, so it should be pro- 
portional to the invariant J = S12 C12 S13 c 2 3 S23 C23 sin 5; 
this invariant can be collected from [C/t\/{7] 3 only. It 
is clear from eqs. (@) and (^|) that Al vanishes in the 
limit r — > 0, so one expects it to be proportional also to 
sin</> = r a /2V a — 2 • 10~ 2 . Therefore, up to a numerical 
constant we have an estimate 

A L -JA 3 sin0. (11) 

The calculations are simplified if |AA/ 2 3 | <C 



\AM{ 



2 i, 1 A A/23 1 ■ ^ n this case, using Eq. 



we obtain 



A, 



[r(j)] 3 

384 



J sine 



h 2 h 2 h 2 . M *2 
AC AB BC 1VI PI 



(AA/ 1 2 3 )i/3(AA/2 



2 2 ) 2 /3 



(12) 



where r(l/3) ~ 2.68 . The estimate ( |l l|) as well as the 
formula ( |l2| ) are valid both in the case when N2 equili- 
brate before tEW and in the opposite case. 

The asymmetry increases when the parameter A ap- 
proaches 1, i.e., the maximal effect for given h a is ex- 
pected when AM 2 ~ (hl b /S) 3 M^f . Notice, however, 
that these and smaller values of AM 2 correspond to sin- 
glet neutrinos strongly degenerate in mass. 

6. Let us present constraints on the parameters of 
singlet neutrinos and discuss the value of the asymmetry. 

(i) To be consistent with the standard mechanism of 
nucleosynthesis, all singlet neutrinos, including the most 
weakly interacting one A3, should decay before the nu- 
cleosynthesis epoch. The decay of iV 3 at T <C Tew oc- 
curs due to its mixing with ordinary neutrino. Requir- 
ing that the decay rate of A 3 exceeds the inverse life- 
time of the Universe at temperatures of order of a few 
MeV, and recalling Eq. (||), we obtain a lower bound 
on the mass, M3 > 1 GeV. Alternatively, for M3 ^> 1 
GeV, we have a lower bound on the Yukawa constant, 



hi > 10- 16 (1 GeV/M 3 ) 3 where h% 



he + S 13^A + S 23^B ■ 
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(ii) If the mass of Na is close to 1 GeV, and/or the 
Yukawa couplings of Na and Nb are close to the bound 
(0), the decays of singlet neutrinos may lead to the re- 
heating of the Universe after the electroweak epoch (but 
before nucleosynthesis), and hence to the dilution of the 
baryon asymmetry. This reheating is rather modest, 
however: given the constraints already imposed, the en- 
tropy density may increase at most by a factor of 10. In 
this case the baryon asymmetry produced before the re- 
heating should be an order of magnitude larger than the 
observed one. 

(iii) Baryon and lepton asymmetries should not be 
washed out before T = T E w by the Majorana mass itself. 
At T ^> Ma the lepton number equilibration rate is sup- 
pressed with respect to the lepton charge conserving rate 
T A , given by Eq. (0), by a factor M\/T 2 , so we have 
to require that tJt E w){M\/T% w ) < H(T EW ), and 
similarly for Nb- The parameters of Majorana singlet 
neutrinos should therefore satisfy 



1 GeV < M a < 100 GeV , 



h A ,h 2 B ^ 1Q- W . 



(13) 



The upper bounds here do not apply to Dirac singlet 
neutrinos. 

The Lagrangian of the model leads, via the see-saw 
mechanism, to the generation of masses of the light (ac- 
tive) neutrinos: m„ a = m a = h 2 v 2 /M a , where v is the 
Higgs vacuum expectation value and v a are mass eigen- 
states - combinations of v e , and v T , The constraints 
@ and (|l3| ) imply that the mass of the heaviest active 
neutrino is in the range m a = (10~ 2 — 10 3 ) eV. From 
the cosmological bound m a < 10 eV we get a constraint 
which is somewhat stronger than Eq. (|l3|), h 2 A < 10 -11 . 
For the lightest active neutrino, the constraints (||) and 
@ lead to m c = (10~ 6 - 10" 1 ) eV. (Notice that the 
mixing parameters of active neutrinos are unconstrained 
by our scenario as they are not directly related to the 
mixing matrix of singlet neutrinos.) 

The above constraints imply that the condition (|K 
is indeed satisfied in large part of the allowed parameter 
space. In particular, it holds in the two cases which we 
now turn to. 

As follows from Eq. (|l2|), in the case when two 
of the singlet neutrinos are relatively strongly interact- 
ing the desired lepton (and baryon) asymmetry ~ 
(a few) ■ 10~ 9 is obtained for generic values of the pa- 
rameters subject to the above constraints. For example, 
for h\~h 2 B = 10~ 12 , AM 2 ~ M 2 and M a = 10 GeV the 
correct asymmetry is generated provided that J > 10~ 3 
which is certainly consistent with Eq. (Q) . The tempera- 
ture of leptogenesis is Tl ~ 10 7 GeV. In this case two ac- 
tive neutrinos are relatively heavy, ttia ~ rns = (a few) 
eV, so that they can constitute the hot dark matter of the 
Universe. Moreover, oscillations between them can solve 
either atmospheric or solar neutrino problem, provided 
that their mass splitting is small. 

A variant of our scenario makes use of weakly interact- 
ing Nq and Nb- In this case two of the usual neutrino 



species have masses in the range (10 6 — 10 v ) eV, and 
the remaining one is relatively heavy. As an example, 
let us take h\ = 5 • 10~ 14 , h% = 10" 15 , h% < h% and 
M a = 20 GeV, which corresponds to the masses of usual 
neutrinos tua ~ 0.1 eV, ms ~ 2-10~ 3 eV and mc < m E - 
This variant fits particularly well into the mass pattern 
suggested by the solar and atmospheric neutrino data 
|L3j. Given that J ^ 10~ 2 due to the constraints anal- 
ogous to Eq. (R), the correct value of asymmetry is ob- 
tained for AM 2 < 10~ 2 GeV 2 . Thus, in this case the 
singlet neutrinos should be degenerate in mass. The tem- 
perature of leptogenesis is lower, Tl ~ 3 • 10 5 GeV. Let 
us note in passing that the degeneracy of masses M a is 
helpful also for obtaining the desired baryon asymmetry 
for very small mixing angles 9ij. 

The crucial feature of the suggested mechanism is that 
it works only if Yukawa couplings of all singlet neutrinos 
are small, h a ~ (10~ 8 — 10~ 6 ). This smallness can be 
explained, e.g., by mixing of N a with very heavy right 
handed neutrinos having Yukawa couplings of the same 
order of magnitude as those of quarks (and charged lep- 
tons) h q . In this case h a ~ hqy/ Ma /Mr, and for our 
values of h a and M a the mass scale Mr may be close to 
the Grand Unification scale, Mr ~ 10 16 GeV. 
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